Surface anchoring on layers of grafted liquid-crystalline chain molecules: A computer 

simulation 
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By Monte Carlo simulations of a soft ellipsoid model for liquid crystals, we study whether a 
layer of grafted liquid-crystalline chain molecules can induce tilt in a nematic fluid. The chains 
are fairly short (four monomers) and made of the same particles as the solvent. They are attached 
to a substrate which favors parallel (planar) alignment. At low grafting densities, the substrate 
dominates and we observe planar alignment. On increasing the grafting density, we find a first order 
transition from planar to tilted alignment. Beyond the transition, the tilt angle with respect to the 
surface normal decreases continuously. The range of accessible anchoring angles is quite large. 

PACS numbers: 61.30.Hn, 61.30.Vx 
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Surfaces orient Hematic liquid crystals, because they 
align nearby particlesau. This phenomenon, called sur- 
face anchoring, isrfipit only interesting from a fundamen- 
tal point of viewoU; jtrfdso plays a key role in liquid 
crystal display deviceaDfi. Surface modifications which 
allow to tune the anchoring properties of a surface - in 
particulafij-the anchoring angle - have great practical im- 
port anceEE!J. 

To be useful in technological applications, such mod- 
ified surfaces must be physically and chemically stable. 
Therefore, preparation techniques where molecules are 
attached chemically to the surface are desirable. The 
coupling between the surface and the nematic fluid must 
be well-defined, i. e., strong. This calls for the use of 
large molecules (polymers), which have themselves liq- 
uid crystalline properties. 

Based on such considerations, Peng, Johannsmann and 
Riihe have recently investigated the anchoring behavior 
of nematic liquid crystals pa .layers of grafted side-chain 
liquid crystalline polymersE^JiliLA " grafting- from" tech- 
nique developed by Riihe et a£§E3 allowed to synthesize 
brushes of up to 230 nm dry thickness. They were forced 
to grow into a homogeneously aligned structure by an ap- 
propriate treatment of the substrate prior to the grafting 
step. When brought in contact with a nematic fluid, 
they transfered their orientation into the fluidli3. Thus 
Peng et al have demonstrated experimentally that sur- 
faces with grafted liquid crystalline chain molecules can 
align nematic liquid crystals. 

The idea of using liquid crystalline polymer brushes 
as alignment layers goes back to Halperin and 
WilliamsLjlij. They suggested to create a competition 
between the orienting effect of the bare substrate, and 
that of stretched chains. The anchoring angle should 
then be determined by the properties of the substrate 
at low grafting densities, and by the chain stretching at 



high grafting densities. Simple mean held considerations 
lead Halperin and Williams to the conclusion that the 
two regimes should be separated by a second order tran- 
sition. This implies that arbitrary anchoring angles can 
be realized at intermediate grafting densities. 

The actual calculation of Halperin and Williams relies 
on two assumptions: The chains are very long and can 
be treated as entropic springs, and they are swollen by a 
good nematic solvent of low molecular weight. At least 
the second condition is not met in real systems: The 
brushes of Peng et al were only swollen to a very limited 
extent, even in a chemically similar solvent. Nevertheless, 
the general idea still seems worth pursuing: Consider a 
substrate which is decorated with liquid-crystalline chain 
molecules (main-chain liquid-crystals for simplicity), and 
which favors parallel alignment: Without grafted chains, 
the substrate aligns the fluid in a planar way. At high 
grafting densities, the chain molecules are forced to stand 
up. The question is: What happens at intermediate den- 
sities? 

The high grafting density limit has been discussed in 
an earlier papercB for the case of short chains in a mod- 
erate solvent. Based on a simple theoretical model, we 
have conjectured that there exists a continuous anchoring 
transition from tilted to perpendicular anchoring, which 
is driven by the grafting density. The theory predicts that 
the anchoring angle (with respect to the surface normal) 
is always smaller than the tilt angle inside of the chain 
layer. Thus the anchoring transition in the bulk pre- 
empts the corresponding tilting transition in the chain 
layer, and may even exist if there is no such tilting tran- 
sition. Computer simulations were consistent with that 
picture, though not conclusive. This anchoring transi- 
tion has nothing to do with the transition predicted by 
Halperin and Williams, which separates a phase with pla- 
nar alignment from a phase with tilted alignment. 
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From a practical point of view, the low grafting den- 
sity regime is more important, since low grafting densi- 
ties can be realized much more easily in experiments. In 
the present paper, we study the anchoring behavior at 
low and intermediate grafting densities by Monte Carlo 
simulations of an idealized, coarse-grained model. The 
general setup of the model is such that, strictly speak- 
ing, the Halperin/ Williams theory does not apply: The 
chains are much too short. Nevertheless, we will see that 
the anchoring transition persists - now as a weakly first 
oder transition. Our results demonstrate that the sce- 
nario predicted by Halperin and Williams is valid in a 
much wider range of parameters than that for which it 
was derived originally. 

In the next section, we introduce the model and de- 
scribe the simulation method. The results are presented 
in secti on We summarize and discuss our findings in 
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II. MODEL AND METHOD 

In our model, both the chain monomers and the solvent 
particles are represented by soft ellipsoidal particles with 
elongation n — cr ond _ cnd /a sidQ _ sidQ = 3. Two particles i and 
j with orientations and Uj (|u| = 1) at the positions 
rj and rj interact with the pair potential 
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where — &o/{ r ij — Oy- + Co) depends on the center- 
center distance r« = liv 



-Tj | and the contact distance aij 



of two ellipsoids in the direction r%j — (r» — rj)/r^-. We 
approximate the latter by the gaussian overlap function 
suggested by Berne and PechukasCJ 
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with x = (k 2 - 1)/(k 2 + !)• The interactions (0) are 
identical for solvent particles and chain monomers. In 
addition, chain monomers are connected by bonds of 
length 6, which are subject to a spring potential with 
an equilibrium length b = bo and a logarithmic cutoff at 
\b-b \ = b s . 
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The chains are stiff by virtue of a stiffness potential 
V A (9i,9aM = -K [ cos(0i) + cos(# 2 ) + 2 cos{9 12 ) ], 

(4) 



which depends on the angles 9± and 9% between the ori- 
entation of a monomer and the adjacent bonds, and the 
angle 0\2 between the two bonds. 

The system is confined between two hard walls at z — 
and z = L z . A hard core potential prevents the particles 
from penetrating the walls, 
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: d z {6) < z < L z - d z (d) 
oo : otherwise 



d z {6) = o-q/2 Vl + cos 2 (<9) (k 2 - 1). 
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The function d z (9) is the contact distance between a sur- 
face and an ellipsoid of diameter ao and elongation k, 
with the long axis oriented at the angle 9 with respect 
to the surface normal. The first bond of each chain is 
attached to one of the surfaces. The grafting points are 
on a regular square lattice. The grafting is fully flexible, 
i. e., no energy contribution couples to the angle between 
the first bond and the surface normal. 

We chose the model parameters k s = IOeo/ctq, k a — 
10eo, bo — 4<7o and b s — 0.8<7o- The simulations were 
performed at the temperature T = 0.5eo/^B and pres- 
sure P = 3eo/<Jo- This corresponds to a state well in 
the nematic phase: The .transition to the isotropic phase 
occurs at the pressureExtEl P = 2.3£o/cq. The bulk num- 
ber density was (p) = 0.313/crg. The length of the chains 
was four bonds, (i. e., four monomers plus the graft- 
ing point) . They are very stiff: The persistence length is 
roughly 210 ao, and the end-to-end distance of free chains 
is 14.8 ao- Compared to experimental systems, our chains 
correspond to main chain oligomers with very few (less 
than ten) mesogenic units. We have not attempted any 
quantitative mapping on a specific system. 

The simulations were conducted in the NPT ensemble, 
at constant temperature and pressure, and a fixed num- 
ber of particles. The systems contained roughly 2000 
solvent particles (the numbers varied slightly in the dif- 
ferent runs) and a varying number of chains with four 
monomers. The simulation boxes were rectangular with 
fixed lateral size L\\ and fluctuating length L z . The lat- 
eral size was kept fixed in order to maintain a constant 
grafting density S. It was chosen Lu = 12uq in most 
runs, except for four runs which sampled grafting den- 
sities between E = 0.111/ctq (4 x 4 chains on an area 
12 • Ylal) and S = 0.174/ct^ (5x5 chains). These were 
conducted in systems with 5x5 chains and slightly larger 
surface areas. The length L z of the boxes was roughly 
L z w 45 — 64ero(±10ero) (higher for higher grafting den- 
sity. Monte Carlo moves were carried out in random 
order and included particle displacements, particle rota- 
tions, and rescaling of the simulation box in the z di- 
rection. They were accepted or rejected according to a 
Metropolis prescription with the effective Hamiltonian 



H = E + PL Z - NT\og(L z ). 



(7) 
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FIG. 1. Total density as a function of the distance z from 
the wall for different grafting densities E. The curves are 
shifted by 0.2/cto f° r better visibility. Grafting densities are 
(from top to bottom) £ a\ = 0.563, 0.444, 0.25, 0.174, 0.111, 
0.0625, 0.0278, 0. 



Here E is the internal energy and N the total number 
of ellipsoids in the system (solvent particles and chain 
monomers). In addition, we have also implemented spe- 
cial Configurational Bias Monte Carlo movescJCj, which 
are specially adapted for our problem: The bonds of a 
randomly chosen chain were removed and redistributed 
such that a new chain was grown from the solvent. We 
calculated the probability P new of constructing this par- 
ticular new chain, and the probability P id of just recon- 
structing the same old chain. The new chain was then 
accepted with the probability 



P, 



min(l,-^iexp(-A£/T)), 



(8) 



where AE is the energy difference between the old and 
the new chain. Details of this algorithm have been pub- 
lished elsewhereo. One "Monte Carlo step" consists on 
average of N attempted translations, 2 N attempted ro- 
tations, one attempt of rescaling the box, and one con- 
figurational bias move. 

Preliminary runs which started from initial configura- 
tions with particles oriented in the z direction showed 
clearly that the walls align the solvent particles in a par- 
allel way. The particles close to the walls turned around 
first, and the parallel orientation then propagated into 
the middle of the film. Unfortunately, this usually lead 
to a state where the director at the two walls pointed into 
two different directions, and was forced to twist around 
in the middle of the cell. In the later runs, we there- 
fore prepared the configurations such that the particles 
all pointed into the x-direction initially. To this end, 
we applied a strong orienting field during roughly 50.000 
Monte Carlo steps. After turning this field off, the sys- 
tem was equilibrated over at least 1 million Monte Carlo 



steps. Data were then collected over 5 million or more 
Monte Carlo steps. 

Due to the complexity of the system, the calculations 
were very time consuming: Every run required several 
months on a Pentium III 800 processor. Therefore, we 
have not yet attempted to study larger systems. Un- 
fortunately, this limits among other the maximum chain 
length - one has to ensure that a chain free "bulk region" 
remains in the middle of the cell. Simulations of systems 
with longer grafted chains will clearly be interesting in 
the future. 



III. RESULTS 

Figure [l] shows profiles of the total density (chain 
monomers and solvent particles) near the wall. They 
exhibit strong layering close to the wall, with decaying 
amplitude further away from the wall. The position of 
the first peak is independent of the grafting density. The 
distances between the other peaks increase with increas- 
ing grafting density. This is a first indication that the 
particle orientation changes from parallel to tilted as a 
function of the grafting density. 

The total excess density at the wall, 

Poxcosb = 77 / dz (p(z) - p bul k), (9) 
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is shown as a function of the grafting density in Figure 
0. Two opposite effects contribute: Particles are pushed 
towards the walls, but cannot move beyond z = <jq/2. At 
low grafting densities, the excluded volume at the wall is 
not quite compensated by the density excess in the den- 
sity peaks close to the wall; at higher grafting densities, 
it is overcompensated. The overall excess density is small 
and rises roughly linearly with the grafting density. 
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FIG. 2. Excess density vs. grafting density E. 
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FIG. 3. Density of solvent particles as a function of the 
distance z from the wall for different grafting densities: (from 
bottom to top) £ erg = 0.563, 0.444, 0.34, 0.25, 0.174, 0.111, 
0.0625, 0.0278, 0.00694, 0. Inset shows a blow-up of the same 
data. 




FIG. 4. Density of chain monomers as a function of the 
distance z from the wall for different grafting densities: (from 
top to bottom) T,crl = 0.563, 0.444, 0.34, 0.25, 0.174, 0.111, 
0.0625, 0.0278, 0.00694. Inset shows the chain density distri- 
butions, i. e., the normalized density profiles p(z)/ J dz' p(z'). 



It is instructive to analyze the profiles of solvent parti- 
cles and chain monomers separately. The profiles of the 
solvent particles (Figure ||) show how the solvent par- 
ticles are successively expelled from the chain region as 
the grafting density increases. The profiles of the chain 
particles are shown in Figure ||. Close to the surface, 
they reproduce the layering of the overall density profile. 
The form of the profiles far from the surface depends on 
the grafting density: At low grafting densities, they de- 
cay smoothly; at higher grafting densities, they develop 
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FIG. 5. Distribution of chain end monomers vs. distance 
z from the surface for different grafting densities: (follow- 
ing the arrows) £ o\ = 0.563, 0.444, 0.34, 0.25, 0.174, 0.111, 
0.0625, 0.0278, 0.00694. Inset shows the same curves on a 
semi-logarithmic scale. 



structure. Additional insight can be gained from looking 
at the normalized chain density distributions at different 
grafting densities £ (Figure [|, inset). They are almost 
identical for the lowest S. As the grafting density in- 
creases, more and more chain monomers are pushed away 
from the surface, i. e., the chains stand up. 

To study this in more detail, we inspect the distribu- 
tions of chain ends in the z direction (Figure |[). At low 
grafting densities, the chain end distributions again re- 
produce the layered structure of the nematic fluid close 
to the wall and decay smoothly, almost exponentially, 
far from the wall. At the grafting density S/ctq ~ 0.15, 
a new broad peak emerges at some distance from the 
wall, which grows and moves away from the wall as the 
grafting density increases further. However, a fraction of 
chain ends still remains in the peaks close to the surface. 
Thus one observes two "types" of chains: Some that lie 
flat on the surface, and some that gradually stand up. 
Increasing the grafting density affects the conformations 
of the chains in two ways: The fraction of chains that 
lie flat decreases, and the chains that stand up extend 
deeper into the bulk. 

Next we study the tilt angle 6 of chains and particles 
with respect to the surface normal. Figure |^ shows the 
average angle (8) of the head-to-tail vector of chains and 
compares it with the corresponding value in a "dry" sys- 
tem of grafted chains without solvent particles. The dry 
layer is a reference system of randomly oriented chains, 
which arc only subject to their mutual repulsive inter- 
actions and to the repulsive interaction with the wall. 
In the presence of solvent, the tilt angle of the chains 
is mostly higher than in the absence of solvent. At low 
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FIG. 6. Average tilt angle (9) of the head-to-end vector 
of the chains as a function of the grafting density E without 
solvent (open symbols) and with solvent (filled symbols). 



grafting densities, the nematic solvent determines the di- 
rector profile close to the surface, and the chains lie down, 
i.e., 9 is large. As the grafting density increases, the tilt 
angle drops rapidly and reaches that of the reference sys- 
tem at S ~ 0.2/cr§. Beyond £ ~ 0.3/cr^, the tilt angle 
almost ceases to decrease. The chains retain tilt even at 
grafting densities where the chains in the dry system are 
almost perpendicular. 

Compared to the dry system, the nematic solvent in- 
troduces a second important feature: It breaks the az- 
imuthal symmetry. In the absence of solvent, the chains 
are disordered. In the presence of solvent, they tilt collec- 
tively in one direction. The direction of the tilt is related 
to the local direction of the director. 
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FIG. 7. Average tilt angle (0) of particles as a function 
of the grafting density E. Open symbols: Chain monomers, 
Filled symbols: solvent particles. 



After having studied the tilt of whole chains, we discuss 
the tilt of individual particles. The average angle (9) 
for chain monomers and solvent particles is shown as a 
function of the grafting density in Figure [?]. The average 
tilt of chain monomers is consistently higher than that 
solvent particles, i. e., the chain region is more tilted 
than the bulk of the system. Looking at the average tilt 
of the solvent particles, one notices a a small jump at 
the grafting density £* w 0.13/<7q. This suggests the 
presence of a first order phase transition. At grafting 
densities lower than £*, the particles lie almost parallel 
to the wall. Beyond £*, they exhibit tilt. 

We can study these phenomena in more detail by 
inspecting the distribution P(cos(9j) of tilt angles for 
solvent particles and chain particles, shown in in Fig- 
ure H It is normalized such that the average number 
of particles in an angle window [9, 9 + d9] is given by 
Ne,de — N total P{cos(9)) sm(6)d9. At low grafting densi- 
ties, the distribution has always its maximum at 9 — tt/2, 
i.e., the particles lie flat at the wall. If one increases 
the grafting density, the distribution first broadens with- 
out shifting. At the grafting density E ~ 0.13/<7q, the 
maximum of the distribution detaches from tt/2 rather 
abruptly and starts moving to lower tilt angles. From 
then on, the maximum tilt angle decreases continuously 
with increasing grafting density. The chain monomer dis- 
tribution has a small secondary maximum at the angle 
tt/2 at all grafting densities. This is the contribution of 
the fraction of chains which are flat on the substrate. 

The maximum angle is plotted as a function of 
the grafting density in Figure ||. The first order phase 
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FIG. 8. Distribution P(cos 6) of tilt angles for solvent par- 
ticles (left) and chain monomers (right) and different grafting 
densities (following the arrows): E ctq = 0.694, 0.563, 0.444, 
0.34, 0.25, 0.174, 0.128, 0.111, 0.0625, 0.0278, 0.00694, and 
(only left figure). The distribution for E erg = 0.128 is shown 
as dashed line. 
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FIG. 9. Position of the maximum angle 8 max vs. graft- 
ing density E for chain particles (open squares) and solvent 
particles (filled circles). 



transition in the nematic film becomes now very evident. 
The state of the film is best characterized by the angle 
distribution of the solvent particles, because the majority 
of them is outside of the chain layer. In contrast, the 
chain monomer distribution characterizes the structure 
close to the surface. The angle 6>max for solvent particles 
jumps to a value different from ir/2 at S* = 0.13/ctq. In 
the case of the chain particles, the 6* max jumps at a slightly 
higher grafting density (X = 0.17/<7q): The chains which 
lie flat strengthen the peak at 6 = n/2. 

Finally, we turn to the discussion of the director pro- 
files. The nematic order parameter and the nematic di- 
rector in a slice of the simulation box are defined as fol- 
lows: From the orientations of all n particles i in the 
slice, one calculates the order tensor, 
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(10) 



where I denotes the unity matrix and (g> the dyadic prod- 
uct. The largest eigenvalue of this matrix is the order 
parameter S in the slice, and the corresponding eigen- 
vector is the director n. The order parameter in a bulk 
system with our parameters is roughly S = 0.75. By 
construction, S increases systematically if evaluated in 
subsystems with few particles. Therefore it tends to ap- 
pear larger in a thin slice than in the bulk. 

In order to calculate order parameter and director pro- 
files, we cut the system into regular slices of thickness 
Az = Icto, starting from the left and the right wall up to 
almost the middle of the film. A thicker slice in the mid- 
dle collected the remaining particles and absorbed the 
fluctuations of the box length L z . This slice was then 
disregarded. 
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FIG. 10. Profiles of the order parameter (S) (thick lines) 
and the components of the director n (thin lines) vs. z in the 
film for different grafting densities E as indicated, thin solid 
line: (|n z |); thin short dashed line: (l^l); thin long dashed 
line: (\n y \). 



Results for four grafting densities are shown in Figure 
|To] . The order parameter S varies very little throughout 
the film. In contrast, the director n changes quite dra- 
matically at higher grafting densities: Close to the walls, 
it is always oriented parallel to the wall, n z -C 1. At low 
grafting densities, it stays parallel everywhere. At high 
grafting density, it turns around very rapidly within the 
first few Co of the surface, and then continues to turn 
more slowly. Substantial reorientation takes place in a 
z-region of about the thickness of the chain layer (cf. 
Figure ^). As a result, the director is tilted in the middle 
of the film. 

The y-component of the director, is overall 

smaller than the a;-component. This is because the 
systems were initially prepared such that the director 
pointed in the x direction. The diffusion time scale for 
reorientation without any driving force is much longer 
than the lengths of our runs (several million Monte Carlo 
steps). In one case, however, we happened to observe the 
onset of a reorientation - the system with grafting den- 
sity X = 0.00694/<7q ventured into a configuration where 
the director twisted between the x and the y direction 
after 5 million Monte Carlo steps. The twist disappeared 
2 million Monte Carlo steps later. It affected the profile 
of (n z ) only slightly. 

Looking at Figure |l^, one notices furthermore that 
the director profile hardly finds enough space to reach 
a plateau in the middle of the film at high grafting den- 
sities. This is a serious problem: Our systems of length 
L z w 50(To(il0<7o) are obviously too small to accommo- 
date a real bulk region between the two surfaces. The re- 
sults are subject to strong finite size effects, and any con- 
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elusions about the anchoring behavior in a semi-infinite 
system must be viewed with caution. Additional finite 
size effects can be expected due to the limited lateral size 
of the simulation box, L|| = 12uq. ft would be desirable 
to study these effects by varying Lu and L z systemati- 
cally in different runs. Unfortunately, an extensive finite 
size analysis is computationally too expensive at the mo- 
ment. Hence our quantitative results can be questioned. 
Nevertheless, we believe that the observed effects will 
persist qualitatively in semi-infinite systems. 

These caveats stated, we collect the results for the 
order parameter n z in the middle of the film at differ- 
ent grafting densities (Figure [ll]) . In a macroscopically 
thick film, it is directly related to the anchoring angle 

via n z = cos (Anchoring). FigUre |ll| confirms and com- 
pletes the picture already suggested by Figures and |: 
At low grafting densities, the surface anchors the direc- 
tor parallel to the wall. The director component (|n z |) 
remains constant and close to zero over a range of graft- 
ing densities, up to £* = 0.13/cr§. Then it jumps to 
a nonzero value and rises continuously thereafter. The 
transition from planar to tilted alignment is clearly dis- 
continuous: A range of values n z is "forbidden" . Hence a 
small range of anchoring angles (Anchoring close to 9 = it/ 2 
(8 e [0.47T : 0.57r]) cannot by accessed at any grafting den- 
sity. As mentioned above, finite size effects have to be 
expected. We cannot exclude the possibility that large 
scale fluctuations drive the transition to second order. 
But since the quantity most susceptible to large scale 
fluctuations in the tilted phase is presumably the az- 
imuthal tilt angle, we suspect that the fluctuations will 
most likely reduce the azimuthal ordering and hence re- 
duce the anchoring angle 9. This implies that the gap 
in Figs. ^ widens, and that the transition becomes more 
strongly first order. 
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FIG. 11. z component of the director in the middle of the 
film as a function of the grafting density. 



IV. SUMMARY AND DISCUSSION 

In sum, we have studied the anchoring of a model ne- 
matic fluid on layers of liquid-crystalline chain molecules, 
grafted on a substrate which favors parallel anchoring. At 
low grafting densities, the substrate dominates and the 
nematic fluid is oriented parallel to the wall. A first order 
transition to a phase with tilted anchoring is encountered 
at a grafting density £*. On increasing the grafting den- 
sity further, the director gradually stands up. The phase 
transition goes along with a major reorganization within 
the chain layer: Below £*, all chains lie flat at the wall. 
Above £*, the chains separate into two distinct classes: 
Some remain flat at the wall, and some stand up. The 
tilt angle 9 of the chains belonging to the latter group de- 
creases with increasing grafting density, and the relative 
number of chains in that group increases. 

The transition is very similar to the anchoring transi- 
tion predicted by Halperin and WilliamsE-3: It mediates 
between a phase with parallel anchoring at low grafting 
density, £ < £*, and a tilted phase at higher graft- 
ing density, £ > £*. It is driven by a competition 
between the anchoring force from the substrate, which 
favors parallel anchoring, and the repulsive forces be- 
tween the grafted chains, which force them to stand up 
at high grafting densities. However, the details of the 
underlying mechanisms are quite different. The theory 
of Halperin and Williams is devised for long chains with 
many hairpins, which can be treated as (anisotropic) en- 
tropic springs. Our chains are too short to support hair- 
pins, and the configurational entropy is much less im- 
portant. The other two constituents of the theory of 
Halperin and Williams, the excluded volume of the chains 
and the elastic energy of the nematic solvent, are obvi- 
ously sufficient to bring about the .phase transition. 

The simulations and theoriesclTLJ consider anchoring 
on layers of main-chain liquid crystalline chains. Exper- 
imentally, layers of grafted side-chain liquid crystalline 
polymers can be synthesized much more easilyBEl. In 
these molecules, the liquid crystalline blocks orient them- 
selves preferably perpendicular to the backbone of the 
chains. In order to generate a conflict between chain 
stretching and the anchoring on the substrate, one needs 
to prepare the substrate such that it aligns the liquid 
crystals in a homeotropic way. So far, this has not been 
done systematically, but preliminary studies by Peng et 
al show that it should be feasible^. The question is: can 
we expect that an effect which has been predicted for 
main-chain liquid crystalline polymers will be observed 
in an experiment with side-chain liquid crystalline poly- 
mers? 

Our study might help to answer this question. Com- 
paring our simulation results with the theory of Halperin 
and Williams, we can conclude that the details of the 
chain conformations and the chain entropy do not really 
matter for the anchoring transition. It is robust enough 
to survive even in a regime for which the theory was never 
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designed. Therefore the chances that it can be observed 
in real systems should be quite high. 
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